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INTRODUCTION 


It Is a well documented fact: that rotating semi-rigid bodies 

("rigid" bodies experiencing energy dissipation) are stable only 

when the rotation is about the axis of maximum moment of inertia. 

This has been demonstrated analytically,^ and. verified in flights 

2 3 

puch as Explorer I and ATS-V satellites . The problem. arises from 

the two potential orientations whii'chT the final spin vector can take 

after large angle reorientation from minor to major axis, i.e., 

along the positive or negative axis of the maximum inertia. 

Reorientation of a satellite initially spinning about the minor axis 

using an energy dissipation device may require that the final spin 

orientation be controlled. Examples of possible applications are the 

4 

Apogee Motor Assembly with Paired Satellites (AMAPS) configuration, 
where proper orientation of the thruster is required; and reorientation 
of ATS-V, where the spin sensitive nature of the despin device ('yo-yo' 
mechanism) requires that the final spin vector point in a specified 
direction. The primary purpose of this work is to investigate 
techniques for eliminating such-spin orientation ambiguities through 
active controls. This control capability would be useful in controlling 
orientation of tumbling bodies for retrieval purposes, and for actively 
reorienting spin stabilized satellites which could be initially spun 
up by the launch vehicle (either earth launch or shuttle deployment) 
about the minor axis of inertia. 
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SUMMARY OF PAST WORK 

Effects of energy dissipation on satellite dynamics have been 

2 3 5 

experienced in actual flights ’ * and discussed in several analytical 
1 6 7 8 9 

treatments. * * ’ * These efforts, have, in general, dealt with 
discussions of the effect of dissipation on the precession rates and 
the total time required for a satellite to go from one spin state to 
another. References 1 and 6 describe attitude drift for single and 
dual-spin symmetric bodies. Reference 7 discusses the effect of 
flexible appendages as energy dissipators on symmetric satellite 
dynamics. Non-symmetric body motion in body fixed and inertial 
coordinates is discussed in Reference 8 while Reference 9 illustrates 
several methods of modeling the energy dissipation process and describes 
effects on non-symmetric bodies. 

Control of attitude drift is discussed in References 10, 11, and 
12. AxisjTnmetric spacecraft control is proposed in Reference 10, in 
which a method is presented for maintaining a satellite on a 
configuration other than spinning about the axis of maximum inertia. 

This is done by addition of energy. into the system to compensate for 
energy dissipation. Reference 11 touches briefly on non-symmetric 
satellite precession control with controlled energy dissipation, and 
presents numerical examples for symmetric satellites. Control of a 
non-symmetric satellite is studied in Reference 12, but only for small 
nutation angles, i.e., active energy dissipation is used to return the 
satellite to the desired stable position after it has been perturbed. 
None of the control systems proposed to date have considered., the 
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question of v/hether the spin vector will align itself with the p.ositive 
or negative axis of maximum moment of inertia when the spin axis is 
initially about the minor axis. 
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DYNAMICS ANALYSIS 

Constant energy motion of a non-symmetric, rigid body, in a 

torque-free environment is geometrically described by Poinsot's motion. 

This is the motion of an inertia ellipsoid (having dimensions 

proportional to the moments of inertia of the body) that rolls without 

slipping on an "invariable plane" while its center is a fixed distance 

9 

/2T/H from the plane. This si t ira tion is illustrated in Figure 1. 

The allowable paths for this motion in the body fixed coordinate frame 
are closed curves called polhodes, each having a fixed discrete energy 
level for given angular momentum. Figure 2 shows several polhodes on 
the inertia ellipsoid for a fixed angular momentum and several different 
energy levels. 

A logical extension of this motion for dissipative bodies is a 

change in the distance from the center of the ellipsoid to the 

invariable plane, resulting in a continuous change in polhode curves. 

Figure 3 illustrates this process along with a geometric interpretation 

of the critical point of the tumbling motion, i.e., whether u has a 

2 

positive or negative component ;,along d^ when the energy passes H / 2 I 2 
determines the final spin orientation. 

Since the nutation angle 9 (instantaneous position of the d^ 
axis with respect to the angular momentum vector) is the parameter of 
concern, equations giving 0 as a function of time and/or rotational 
energy are desirable. Although the equations of motion for a non-sym- 
metric body are elliptical in nature and, therefore, not readily 
reduced to an analytic solution, limits on the nutation angle can be 
developed as straightforward functions of the energy state. For 
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Figure 1. Rigid Body Poin'sot Motion 
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Figure 2. Rigid Body Polhodes 
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Figure 3. Dissipative Body Polhode 


8 


T < n Reference 9 develops these bounds in the following manner. 

Extreme values of 0 exist when OJ lies in either the 2-3 or 1-3 plane, 
as explained in Appendix A and illustrated in Figure A. This results 
in -the angular momentum and energy taking on the following values : 


at 6 


at 6, 


„2 ^2 2^^2 2 

H +l3<.3. 


„2 ,22^,22 

H - I3 0J3 + I3 «3, 


2 T . 130 , 3 ^ + 13.032 


2T = 130,32 ^ 3^^^2 


From the geometry of the situation, it is evident that 


T 2 2 

, 2 . h “2 

sin 0 == z — 


T 2 2 

2. ^1 “1 

sin 0. = z — 


H 


I (21 - I 3 OI 32 ) 


i <22 - 

n 


I 3 H 


2 222 

T < 213 ^ - h" + 




1 2 2 2 
^ ( 2 I 3 T - ) 


Manipulation of these equations ‘leads to 


- 1,(21 T-H^) 

sin^e = -j 3 ^ ^ 

“ H (I 3 -I 2 ) I 3 H (I 3 -I 2 ) 


{ - 2I3T + CI3-I2) ^ 


. 2. Ii(2I3T-r2) 1^2 ^ ^ ^ ^ 

sxn 0- - —z + Y 


ji— ( - 2I3T + h2 + I30,3'(l3-l3) ) 
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Figure 4^ Inertia Ellipsoid Orientation at Extreme Values of 0 
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Keeping in mind the particular energy and momentum conditions that 

exist at the respective extremes, the terms in brackets reduce to 

zero, leaving . 

. I (21 T-H^) . 

ain^e = -4 


2 I (21 T-H^) 

= i— ^ — 


a‘(i3-ii) 


With the energy above the separatrix level (T > H tll^ the following 

parallel development can be made. Figure 5 shows the inertia 

ellipsoid with the d^ axis in one of its extreme positions. As 

2 

never equals zero (the elliptic solution for 03^ with T > H never 

13 

changes sign ), the upper and lower limits must both exist where 0)^ 

equals zero. At this point o)^ and determine the nutation angle 0. 

For a given energy state, every zero of the function coincides with 

the same value for with having the same absolute value but 

alternating sign. Thus the d^ axis must oscillate about 90° with the 

maximum angle of the oscillation giving both the upper and lower 

bounds on 9. Defining this angle., as 0 results in the following 

ex 

equations : 


2 2 

. 2 „ ^3 ^^3 ^3 , 2 . 

®ex = — (2T - I w ) 

® H H ^ ^ 


IlH 


^ (2Ij^T - + 1^0)^^) 



Inertia 

Ellipsoid 



Figure 5. Extretne Position of d Axis 
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The above relations take into account the fact that at this point in 
the motion 

2T - 


„2 ^2 2^^2 2 
H = 1 ^ 03 ^ +13033 


Further manipulation yields 


^2 = 0 


2 I (21 T-r) 

sin 0 = -4 + 




^ { 13^3^(13^-13) - (2I^T-H^) > 


Again, bearing in mind the particular energy and momentum conditions 
at this point, the final term in brackets equals zero. Therefore, 


2 21 I I 

sin‘^0 = ^ — T 

ex 




To provide as general an interpretation as possible, introduce the 
following terms; 


h2 - 17 


T* = 


max 


where T^_^ is the maximum possible energy state for a given angular 

2 

momentum. Thus, T = H /2I. . Making the above substitutions, the 

max 1 

bounds on the nutation angle reduce to the following forms : 
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a. Motion above the separatrix (T > H /2I^) 


0 = 90° + arcsin 

ua 


1-T 


13 


=90 - arcsin 


/ l-T* ‘ 

/ 1-Ii3 


b. 


Motion below the separatrix (T < H /ll^) 


9 , = arcsin 
ub 


T*-I 


13 


’■12 ’•13 


'T*-I 


= arcsin 


13 


l-I 


13 


At this stage, two relationships are to be noted: 

1) The intermediate moment of inertia (contained in the term 
^12^ has no effect on the nutation angle range above the 
separatrix, and no effect on the lower limit of 0 below the 

t 

separatrix, 

2) At the separatrix 

T 


Ji 

21 


I* = — = I 

I2 12 


For a specific body having a minimum to maximum inertia ratio of one 
half ” 0*5) the range of T* is given by 0.5 < T* < 1.0. Figure 6 

shows a plot of the bounds on 0 vs T*, assuming that " ^13' 
.Therefore, the first pair of equations holds for the entire range 
of T*.. This case corresponds to a thin rod, initially spinning about 
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its axis of minimum moment of inertia. As any transverse axis has 
the same (maximum) moment of inertia, the final nutation angle can 
take on any value between 0° arid 180°. Figure 7 shows the bounds on 
6 for a body having the same 1 ^ 3 * ^ intermediate ratio, 

1^2 “ 0.8. This value for T* (T* “ lj^2 ~ 0*8) can be seen as the 
critical point in the motion, as the bounds on 9 divide into two 
regions: one ending in a final nutation, angle of zerp;, the other 

ending in a final 9 of 180°, i.e-;* d^ pointing in the opposite 
direction as (o and h. 

Noting the continuity of the minimum values for 0 at the 
transition point, the following relations can be written: 


6, - 90° - arcsin 

la 1-1^3 


= arccos 


1-T* 


l-I 


13 


T*-I 




13 


l-I 


13 


T*-T 

2_ 1-T* ' . 2. 13 

cos 9„ = • • V ' ,"'sin 0„, = -r-z 

Aa 1~^13 


2 2 ^”^13 

cos 0- + sin 9. = r— — 

x,a Jib ^“^13 


= 1 


®Jia ®Jlb 



iso'f 
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(31in.T:'i jaMOT puE jaddf]) 0 
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This Is consistent with the earlier observation that 1^2 

affect the bounds 0„, , 0 and 0. , as the value for has no 

£b ub £a 12 

effect on these portions of the curve. 

From the above description of tumbling motion it is apparent 
that the critical point in the motion is where the energy satisfies 

T = H^/2l2 . 

At this point either the nutation angle or the sign of identifies 
which direction the final spin vector will assume. As the nutation 
angle is rather difficult to measure, the sign of ti)^ will be used to 
determine the need for a control device. Assuming the control 
algorithm outlined in Figure 8, magnitudes of the three components of 
6), along with the sign of must be known. Sensors are readily 
available for such application. Thus, it will be assumed that such 
equipment will be included in the satellite. 

With the rotation vector o) known, real time calculations of the 
energy state and the dimensionless energy parameter T* can be made. 
Comparison of T* with the constant 1^2 then determines when the 
critical energy state has been reached, as T* = 1^2 ®t the critical 
point. 
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Figure 8. Algorithm for Elimination of Spin Ambiguity. 
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CONTROL SYSTEM ANALYSIS 


Assuming that some control is required, three methods are 
discussed here. Moving masses in such a way as to lower the ratio 
1^2 results in .T* again being greater than the naotion 

continues with oi lying above the separatrix level. When becomes 
positive, the masses are returned to their original positions, and 
1^2 is raised to its original value. While some energy is added in 
this process the net effect still leaves u) crossing the separatrix 
in the desired area. Starting an electric motor within the satellite 
increases the energy, and therefore T*, again allcn^ing more time for 
0 ) to move around until is positive. At this point, the motor is 

stopped, and again w crosses the separatrix in the desired area. 


Both of the above systems are attractive in that they are completely 


internal in nature. Thus, the angular momentum vector remains 
unchanged, i.e., given the initial conditions, the final state is 
dictated. However, the size and orientation of the control elements 
must be such that sufficient change in 1^2 "T* roust be generated 

to insure that w does not cross the separatrix a second time before • 
0)^ becomes positive. Depending on the size of the satellite and the 
energy dissipation rates, this may require excessive size or volume 
requirements to effectively control the spin. 

In these cases it may be necessary to go to an external torque 
to control the spin, i.e., reaction jets. These could be applied to 
increase the energy should the need arise, and decrease it again at 
the appropriate time. This, however, will result in some shifting of 
the final angular momentum vector. A final choice of a control system 
will depend on the individual situation. 





20 


For the first method suggested, motion of the masses could 
result in some loss of energy (see Appendix B) . Thus, some criteria 
are required to verify that the change in is sufficient to over- 
come any loss in energy. For the configuration shown in Figure 9, 

=1, + A ■+ 2mr^ 

J.XJ X 

= l^ + -& + 2mR^ 

13 ^ = 13 + 0 + 2 m (R^ + r^) 

where the subscript *B' denotes before the weight shift. Similarly, 
with 'A* denoting after the weight shift, 


^lA - 


IjA • ^26 


Sa ' Sb + 


From these, the following momentum and energy relationships can be 
written: 

\ " ‘‘^IB^IB^ ,‘^2B^2B^ “3B^3B^ 

17 7 2 

T = — {l w + I w + I 0 ) } 

B 2 '• IB IB ^25 2B 3B 3B 


*B 


maxB 21 


IB 


•j* = 


B T 


maxB h 


iS-fT 2^_ 2^^ 2, 

2 ^^1b“iB ^2b‘^2B ^3b''^3B ^ 

B 
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2A 2A 


3A 3A 


2 ^^1a‘^1A ■‘’ ^2a“2A "'■. ^3a“3A ^ 


maxA 21 


lA 


"*A ■ f W ^ ^2a“2a' ^ ^3a“3a' > 

A 

The magnitudes of h. and h must be equal. If the masses are assumed 

A B 

to move quickly xfith respect to the rotation rates of the satellite, 
the components of h. can be assumed to be equal to the components of 

A 


hg. These conditions yield: 


•“a " *‘b * ’> 


IB 

W, . = rz — 03, „ 
lA I, . IB 
lA 


12B I 


^IB 

2B 


2B 

“2A “ l2A 


3B 

0 ) = 03 

3A I_. 3B 
3A 


12A I 


14 

2A 


At the crossover point, T* = 

B xZB 


T 2 ^ 2 . ^ 2 h 

^1b‘*’iB ^2b“2B ^3b“3B I 


2B 
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The desired result is for T*. > I,-. 

A IZA 


^lA 2 ^^ 2 _^^ 2 ,- ^ 

.2 ^^lA^^lA ^ ^2a“ 2A ■*■ ^3A 3A ^ I- 
h , 2A 


2 ^ ^ 2 ^ 2 ' h 

^lA^^lA ^2A^^2A ^3 a“3A I,. 

2A 

2 

Substitution of h from the equation for. T* at the crossover point 

B 

results in: 

2 r^lB^2A“’’2B^lA _ 2 r^3B’'2A"^2B’'3A , . » 

‘“u '^2 'a ■■■ ' " ^3b“3B 

To insure that this is true (thus T*^ > I]^2A^ * 


^1B^2A ■ ’■2B^1A ^ ° 


^3B^2A " ’•2B^3A ^ ° 


Substituting for I 2 A* ^3A 


I^b( 2R + r) + r I^j > 0 


IjB(2R+ r) -Rl2B>0 


The first condition is always satisfied. The second can be 


2+i>^ 

^3B 


rewritten as: 
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As I„n < I-D (by definition of the problem), this condition is 
satisfied for any positive values for r and R. Thus, the motion of 
any masses, in the configuration illustrated in Figure 9 will result 
In the motion subsequent to the weight shift being above the 
separatrix (T*^ > ' 

Use of the momentum wheel, the second method proposed, is the 
simple addition of energy to the system.. For the- simulation case used, 
the motor was aligned with the axis, to produce as large an effect 
as possible on the angular velocity vector to. The reasoning behind 
this choice is shown in Figure 10 where it can be seen that the 
easiest way to physically move to above the separatrix is to increase 
(Oj^, Thus, the control motor for the simulation was mounted 
accordingly. 

The final method considered is the straight forward application 
of an external torque through the use of attitude control thrusters. 

As mentioned previously, this is probably the least desirable of the 
potential control methods, as it results in a shifting of the angular 
momentum vector in inertial space. However, it is a possible means 
of controlling the final spin by the simple addition of energy 
(positive torque) should OJ cross the separatrix at the incorrect 
interval, and the reduction of energy (negative torque) when the 
desired state has been reached. 
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SI^^^JLAXIOI'!S AND KE3ULTS 


The dynamic simulation used is very close tp that used in 

Reference 4. This consists of a rigid satellite with a fluid ring 

d^per mounted around the axis of maximum moment of inertia, in the 
* 

plane of the transverse axes of inertia. Partially filling this ring 
with damping fluid results in some motion of the center of gravity of 
the satellite. Use of the characteristics for the partially filled 
damper with Euler's equations (for motion about a fixed point) results 
in some oscillation in the inertial angular momentum. Since the 
principal concern of this study is the relationship between the 
spacecraft orientation and the inertial angular momentum, two alterna- 
tives were considered for correcting this problem: (1) use of 

equations derived in Reference 14 for motion of a body with moving 
masses, and, (2) use of a fully filled fluid damper, eliminating the 
moving mass center. Use of the first method would involve considerably 
more computer time to handle the complex differential equations arising 
from the motion of the center of mass. This added time and effort 
would add little or no information concerning the point of interest 
here, control of the final spin vector. The se'cond method, while 
seldom used in practical applications due to very low dissipation 
rates as steady spin is approached, does provide realistic energy 
dissipation in the region of interest, with the advantage of a simple, 
relatively accurate mathematical model. Thus, the second method was 
chosen and the equations of motion derived accordingly (See Appendix 
C) . For the weight shift control method, these equations (before the 
weight shift) are; 
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0)1 = 


102^3(^2 ^ ~ ^3 “ ^ “ Co) 2 (W 3 + a) 

II + A + 2mr^ 


0)„ = 


0)l0)^(ll + 2mR - Ii - A) + Co)i(o)i + ct) 
I 2 + A + 2 mR" 


- ^ 0.C395 CR jP 1 • I o 

2 ^2, , , . ^ 2 


a)iW2(Ij^+2ni(r -R )-l2> + 


SR a 

n> 




‘^3- 2m (R^ + r^) 


.. -0.0395 R ,P i;i - 

ad ^ a U * 2 

m — ^ " “3 

SR ■' a ^ 


Appropriate changes in the moments of inertia and application of the 
same manipulations provide similar equations for the motion after the 
masses have been shifted. These equations are then solved in the 
simulation diagrammed in Figure 11. For all the simulations, a 
fourth order Runge-Kutta method is used to solve the differential 
equations, using a time step of 0.25 sec. Several time steps were 
tried, but 0.25 sec. gives the best compromise between desirable 
accuracy and reasonable simulaSison time for a given amount of computer 
time. 

All equations are solved in the body fixed coordinate system. 

To provide the nutation angle 0 and to verify the accuracy of the 
entire simulation, Euler's angles ((J»,<}),6) are calculated using the 
following equations. 


■ 4 ^ = 77 


sin 0 ' ^“1 “2 



27 



Figure 11. Flow Chart for Computer Simulation 
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• COS 0 

4^ *= oj- - — r— H" sin <}) + w cos ({)) 

j sin-. 0 1 ' i. 

0 = cos 0 - (jO^ sin (j) 

again with a fourth order Runge Kutta -method. The values for Euler’s 
angles are then used to transform the. angular momentum from the body 
fixed to inertial co-ordinate systems. As this should remain 
constant -for .the first ..two .control -systems, it provides an easy check 
on the validity of the equations and the adequacy of the time step. 

The spacecraft modeled is the vehicle presented in Reference 4. 
The appropriate physical properties are; 

= 1785 ft-lb-sec^ (2380 kg - m^ 

I 2 = 5190 ft-lb-sec^ (6940 kg - m^) 

= 6920 ft- lb-sec^ (9250 kg - m^) 

For the fluid damper the following figures apply: 


p = 1.27 X 10 ”^ ■ — y- - (1.36 X 10^ kg/m^) 

In^ 


V = 1.8 X 10 ^ in^/sec (1.16 x 10 ^m^/sec) 


R , = 4.35 ft. (1.33 m) 
ad 


mercury as 
viscous fluid 


D = 0.13 ft. (0.0398 m) 

For the control masses, the following values were chosen: 


R = 1 ft. (0.30 m) 


r = 3 ft. (0.91 m) 
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m = 5 lb. (2.26 kg) 

It should be noted that simulations were run with smaller lengths 
and lighter masses. However, these proved insufficient for having 
the desired effect. With the above figures, the parameters 1^2 
take on the values 0.402 and 0.289, respectively. Figure 12 
shows the theoretical limits on the bounds on 0 for this particular 
body. Figure 13 shows the results of the computer simulation for 
the angle 0. This case considers an angular momentum equivalent to an 
Initial spin rate of ~ 10 rpm. The initial conditions are chosen 
such that the satellites energy state is relatively near the 
separatrix, yet far enough away to allow observation of several 
cycles before the critical point is reached. It should’ also be noted 
that initially the fluid slug has a body relative velocity of zero. 
Figures 13 and 14 show the nutation angle, energy state, and 
dissipation rate for the first simulation. As can be seen from 
these figures, no control is required for this case, as 0 is oscillating 
tCT/ards the final limit of 0°, and is oscillating towards a final 
positive value. These calculations are based on the empirical 
friction coefficient f = 0.316/R for the fluid damper. Changing 
this coefficient by 5% has the effect of reversing the final outcome 
(see Figures 15 and 16). Thus, any attempt to predict the crossover 
point and control the final spin orientation by alteration of the 
initial conditions would require knowledge of the energy dissipation 
rate, to within at least 5%. Alternatively, Figure 17 shows this 
sane second case with the moving masses incorporated into the system. 

As the energy reaches the separatrix level initially, 0 is greater 
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rigure 13. Nutation Angle and Energy for Desired Result 
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Figure 14. Dissipation Rate and a)_ for Desired Result. 


Simulation results for 



Figure 15. Nutation Angle arid Energy for Incorrect Result 



Figure 16. Dissipation Rate and for Incorrect Result 






Critical Point-— V > ^Masses returned to 

(^aasses moved) \ original positions 



Figure 17. Nutation Angle and Energy Parameter for Moving Mass Control Method 
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than 90°, therefore, it would oscillate towards the final va^lue of 
180°. Movement of the masses at this point results in the lowering 
of the ratio therefore the critical energy level, allowing 

0 to oscillate below 90°. At this point, the masses are returned to 
their original positions returns to its original value. The 
energy added by the movement of the masses is quickly dissipated, and 
the critical energy is again reached, with 0 now oscillating towards 
0 °. 

Using the electric motor control method and the initial conditions 

used to obtain Figures. 15 and 16, the simulation produced the results 

shown in Figure 18. Rotor characteristics were taken as 

=0.05 ft-lb-sec^ (0.067 kg - m^) 

*2R “ ’•3R ” ft-lb-sec^ (0.034 kg 

h„ = 7.5 ft-lb-sec (10.0 ) 

R sec 

these being the characteristics of a large, commercially available 
momentum wheel. Three such wheels were required to produce the 
desired results. Figure 18 shows the energy state of the satellite, 
with addition of the rotor energy an obvious platform in the energy 
level as the critical state is reached with 0 > 90°. Shutting off 
the motors as 0 reaches 90° again allows 0 to oscillate toward a 
final value of 0°. 

Using attitude control thrusters, assuming a thrust of 5 lbs. 
for 3 seconds at a radius of 6 ft. produces the results shown in 
Figure 19. Again, the activation of the control system results in an 
increase in the energy, allowing 0 to again return to the desired 
region (0 < 90°) . Theoretical limits on 0 are not given for this 
case after the application of the control torque, as the theory gives 



Initial Critical Point 



Figure 18. Nutation Angle and Energy State for Electric Motor Control Method 
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0 as the angle between and the angular momentum vector. In the 
preceding cases this coincided -with the z-axis of inertial space. 

The external torque results in a deviation of the direction of the 
angular momentum vector of approximately 3.4° while 0 is still being 
measured from the z-axis. This accounts for the irregular nature of 
the oscillations after the application of the control torques. 
Application of the external torque also results in a change in the 
magnitude of the angular momentum vector of approximately 1%. 
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CONCLUSIONS 


The analytical expressions developed for the nutation angle 0, 

coupled with the simulation results coinciding with the predictions, 

appear to verify that a tumbling body exhibits two distinct types of 

2 

motion. For an energy state greater than h / 2 I 2 (where ^ 3 ^ ^2 ^ ^3^ 

the nutation angle oscillates symmetrically about 90°, while 

oscillates through positive and negative values. At the point where 

2 

the energy equals h / 2 I 2 motion is transformed, with 0 and still 

oscillatory functions; however 0 now oscillates toward either 0° or 

180°, while remains either a positive or negative function, 

oscillating towards the final stable value in either the positive or 

negative sense. Control over this final orientation can be most 

2 

easily accomplished at the critical point T = h by varying 

either the moments of inertia, the energy state, or both. The work 
presented here verifies that the motion can be controlled by the use 
of either moving masses, electric motors, or thrusters, along with 
minimal sensing equipment. The most attractive method appears to be 
the moveable masses, as relatively large electric motors were required 
to handle the sample case presented. Thrusters, while originally felt 
undesirable, may, with further work, prove to have one large advantage 
This lies in flexibility of control. While either the mass or the 
motors would have to be sized for worst case conditions, it may be 
possible to derive an optimum thrust for any particular case, thus, 
minimizing both the deviation in the angular momentum vector and the 
mass needed to effectively control the motion. 
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Figure A-1 
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Extreme values of 0 coincide to points where either 


equals zero 
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APPENDIX B 

Loss of energy through mass motion 



Consider a simply rotating body as 

Illustrated with moment of inertia I. 

The two masses at distance fi. from center 

2 

result in a total = I + 2mi 

liierefore, the momentum and energy equations are 

h = w(I + 2m!L^) 

T = Y (I + 2ml^)i/ 

Allowing the masses to move out to a length 25., the moment of inertia, 
momentum, and energy become: 

= I + 8m5,^ 
h = 0)^ (I + 8m5,^) 

Tj = j (I + Sms}) 

where = final rotation vector. To maintain angular momentum, 

,I + 2m5.^. 

U, = ( r) W 

I + 8m5, 


_ 1 ,I + 2m5. ^ „ „2^ 2 

T, = 7 ( r) (I + 2m5, ) (0 

I + 8m5, 


I + 8mZ 


Tj < T. 



of 


APPEiVDix C 


equations /, 

of motion^ 


"i ' “l»l + A + 

"2 ' “2«2 A + 2^2, 

"3 ■ “3«3 + 2„ f,2 , 

+ CC'-j + i) 

“i ““i«j+A + 2 „ 2 ^ 
'’2'“2«2*A+2«2^ 


Euler* 


'"a “ + 


equations i 5 




yield 

+ A + ,„,2, 


* + cci + 


!3 + a) 


1 ^ + 2 mr') + ,, 

“ 2"3 a, + 


3 ”3 + 2 a (^2\ 2 


* ^3 CSo.fu 


' 2 '“3 + CO 


■ “^“3 «2 - A . 2a.3, . „ 


“2<'I, + A + 2„o2 


2 oiR^) . 

■' ^ (t},U) fT . 

i 3 + A + 


- »,<«. 


i“3 «3 3-Cr^ 


f "3 ^c:) =0 


“3^^3 + 2m(r- 


■'■ + r/' 


i“2«2 + A + 2 „r2j 


1 + A + 2n,r2) , ^ 
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+ A + 2rar^) = ^ “ 2mr^} - Cni^ (w^ + “) 

\0-2 + A + 2mR^) = {l^ + 2mR^ - I^ - A} + Co)^ (w^ + a) 

w^Cls + C + 2m(r^ + R^)) = {l^ + 2m (r^ - R^) - I^) - C a 


To find an expression for a, consider the friction force on the wall 
as in elementary pipe flow. « 

« (f/4) pv^/2 


where *f is an empirical resistance coefficient. The torque then 
becomes: 


N. 


- *ad ( 


R 


a 


where the term in parenthesis insures that the torque acts opposed 
to fluid slug motion. Using an empirical friction coefficient by 
Blasslus of 

f = 0.316/R 
n 


the torque then becomes: 


Letting 


-0.316 pv R ,P3 , , 

N ^ 

3 1/4 '• • 

8R a 


V = R, a 

ad 

R = R . a D/v 

n ad 
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-0.0395 p 

a (“) 


(D/V) 


IM 


N, 


-0.0395 p 


Writing the kinetic energy of the fluid slug as 

1 2 1 2 1 • 

Tg = "2 + 2 ^2 ^ 2 ^ (w^ + ct)" 

a Lagrangian formulation can be used to relate a to the motion of the 
body. In the following equation, is the retarding torque, therefore 


‘!a-''3 


i_ /M) . 51 , Q 


da 


3a ^a 


.0395 |a| 


(C(u )3 + a)= - 


1 11/4 


-.0395 p 

CO - 1/4 - CU3 


(D/V) 


a 


a 


-.0395 p 

C(D/v)^^^ a 


- 


\ 
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Gomb'inlng this with the equations for u results in the complete 
equations of motion for the rigid body 




{I2 ^ “ ^3 “ ^ " 2mr } - Cui^ct 

+ A + 2nr^ 


“2 “ 


(OlWa {I3 + C - I3 - A + 2mR } + Ciii^a 
I2 + A + 2inR^ 

- , .0395 pR ,P3 ]a| 

0) w. {I + 2m (r - R^) - I } + - \ 7T \ 

I, + 2m(r" + R^) 


11/4 


.0395 p R .P3 

.. ad ' ‘ • 

a = TT7~. ‘^•5 

C(D/v)^^^ a ^ 

• 

For the computer simulation a is created as the principle variable. 
Thus the preceding are four first order differential equations in 
four unknowns. Moving the masses to the ^ axis (the proposed control 

motion) jresults in the following equations; 

^ W2W3(l2 + A - - C) - Cco^ a 


+ A 


+ C - I3 - A + 2m(r + R)*^) + a 


“2 ” 


I2 + B + 2m(r + R)‘ 
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